Expansion of a quantum degenerate boson-fermion mixture 
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We study the expansion of an ultracold boson-fermion mixture released from an elongated magnetic 
trap, in the context of a recent experiment at LENS (G. Roati et al, Phys. Rev. Lett. 89, 150403 
(2002)). We discuss in some details the role of the boson-fermion interaction on the evolution of 
the radial-to-axial aspect ratio of the condensate, and show that it depends crucially on the relative 
dynamics of the condensate and degenerate Fermi gas in the radial direction, which is characterized 
by the ratio between the trapping frequency for fermions and for bosons. Our numerical simulations 
are in reasonable agreement with the experiment. 
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I. INTRODUCTION 

The experimental realization of ultracold bose-fermi 
mixtures of alkali atoms introduces a novel tool for the 
study of various quantum phenomena [1-4]. The most 
appealing one is certainly the Bardeen-Cooper-Schrieffcr 
(BCS) type superfluidity, since the presence of boson- 
fermion interaction in the mixture is expected to in- 
duce an effective attraction between fermions by exchang- 
ing density fluctuations of the bosonic background [5,6], 
analogously to that of superfluid 3 He and 4 He in low- 
temperature physics [7,8]. With this respect, the recently 
reported quantum degenerate mixtures of 40 K (fermion) 
and 87 Rb (boson) by the LENS group would be par- 
ticularly interesting due to the large boson-fermion at- 
traction, which has been directly manifested by the col- 
lapse of the degenerate fermions above a critical number 
of particles [4] , or alternatively, by the observation that 
the Bose-Einstein condensate (BEC) coexisting with the 
Fermi gas inverts its aspect ratio more rapidly than a 
pure condensate during the ballistic expansion [3]. In 
the latter case, the experimental result agrees qualita- 
tively with the expectations of a tighter confinement for 
a BEC in a Fermi gas with mutual attraction [9,10], and 
is qualitatively fitted by the theoretical predictions for 
the free expansion of a pure BEC with trap frequencies 
10% larger than the actual one. However, possible effects 
of the boson-fermion attraction during the early stages of 
the expansion are not considered in the fitting and cer- 
tainly need a further investigation. 

In this paper, we would like to study in details the 
effect of the boson-fermion attraction on the expansion 
of the condensate as well as of the Fermi gas after the 
release from the strongly elongated magnetic trap. We 
have taken into account both ground state effects and the 
effect of attraction between bosons and fermions during 
the early stages of the expansion. The latter one is found 
to play an important role on the evolution of the radial- 
to-axial aspect ratio of the condensate. Our numerical 



simulations are in reasonable agreement with the exper- 
iment. 



II. FORMULATION 

We consider a dilute spin-polarized boson-fermion bi- 
nary mixture at very low temperature trapped in a 
strongly elongated harmonic oscillator potential. In the 
semi-classical Thomas-Fermi approximation, the conden- 
sate and degenerate Fermi gas evolve according to the 
Stringari's hydrodynamic formulation [11] 

— — + V • (n b v h ) = -T coU , 
at 



m b - 



dt 
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and Boltzmann-Vlasov kinetic equation [12,13], 
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respectively. Here /(r,V/,t) is the single particle 
phase space distribution function for fermions; n b (r,t), 
iif (r, t) — j d 3 Vff (r, v/, t) and v& (r, t), v/ (r, t) are the 
densities and velocity fields for bosons and fermions, re- 
spectively. Here, for simplicity, we consider the case in 
which they are trapped by symmetric confining poten- 
tials, V%? (r) = \m b j (t^xbjP 2 + ^6,/ z2 )' in a concen- 
tric configuration. Note that in the realistic experimental 
situation the centers of mass of the condensate and Fermi 
gas are instead displaced due to the different gravitational 
sag for 40 K and 87 Rb [3]. The displacement is not suffi- 
ciently large to affect the geometrical overlap of the two 
degenerate species. However, in some degree it reduces 
the boson-fermion attraction. Of course, a quantitative 
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analysis should take into account this displacement. We 
will comment on this point later. 

The influence of the interatomic scattering processes 
between the two species on their dynamics has been 
taken into account in two parts: mean fields and col- 
lisions. In the above equations, the term gbfTif (r, t) 
and gbfnb(r,t) are included as the Hartree-Fock mean 
field effect of boson-fermion interaction [14] (note that 
the latter one is known as the Vlasov contribution in 
the literatures) . The boson-boson and boson-fermion in- 
teraction strength for the pseudopotentials, gbb and gbf, 
are related to the s-wave scattering lengths abb and abf 
through gbb = ^h 2 abb /nib, gbf = 2irh 2 abf /rribf, in which 
mbf = mbmf I (nib + m f) is the reduced boson-fermion 
mass. On the other hand, the dissipation term — Y co u or 
I co ii accounts for collisions [15]. The relative importance 
of the mean field and dissipation terms in describing the 
dynamics of whole system depends upon the rate of col- 
lisions compared to the characteristic time scale of the 
single particle oscillation. When collisions are relatively 
rare, the mean-field interaction dominates and the sys- 
tem is in the collisionless regime. In the opposite, hy- 
drodynamic regime, the degenerate Fermi gas is always 
in the local equilibrium with the condensate [16]. The 
dynamics is then most appropriately described by the 
hydrodynamic Euler equation of motion [17,18]. Due to 
the lack of information on the collision rate in experi- 
ment, we shall mainly focus on the collisionless regime 
by neglecting the dissipation terms — T co u and I co a . The 
opposite hydrodynamic regime will be briefly discussed 
at the end of this paper. 

In the absence of the boson-fermion interaction (gb / = 
0) and dissipation terms, both the Eqs. (1) and (2) admit 
a simple scaling solution, i.e., 



n b (r,t) 
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for the condensate [18,19] and 



(3) 



/( r ,v/,i) = /o 



7i(*) 



,V(r,t) 



V i (r,t)= li (t)v fi -^r i 



for the degenerate Fermi gas [12,13,20]. Here n° and /o 
are the equilibrium distributions. The dependence on 
time t is entirely contained the six dimensionless scaling 
parameters, hi (t) and 7i(t), where i = x, y, and z. By 
substituting this solution into Eqs. (1) and (2), it is easily 
to show that, the scaling parameters obey the differential 
equations [19,20], 
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Solutions of Eqs. (5) and (6), respectively, determine the 
time evolution of the pure condensate and Fermi gas. In 
particular, they can be used to study the expansion of the 
system after a sudden and total opening of the trap at 
t = 0. For an elongated cylindrical trap with anisotropic 
parameter A = oj z b/uj±b (= u z f/u±f) <C 1, one may find 
that [19,20] 



6±(t) = Vl + r 2 , 

b z (t) = 1 + X 2 t arctan r — In \J\ + t 2 , 



and 



7±(t) = Vl + ^r 2 , 
7»(r) = 



(7) 
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where we have introduced a dimensionless time variable 
t = (J±bt and (3 = u±f/u)±b- 

In the presence of the boson-fermion interaction (g^f ^ 
0), however, the simple scaling solution is no longer sat- 
isfied at every position r after the substitution. A useful 
approximation, in the first order of gbf , is to assume the 
scaling form of the solution as a priori, and fulfill it on 
average by integrating over the spatial coordinates. This 
strategy has been recently used by Guery-Odelin [12] to 
investigate the effect of the interaction on the collective 
oscillation of a classical gas in the collisionless regime 
and by Menotti et al. to study the expansion of an in- 
teracting Fermi gas [13]. In some sense, this approxima- 
tion is equivalent to the sum-rule approach [13] that is 
extensively used in evaluating the low-energy collective 
modes of dilute quantum gases [11]. We have also re- 
cently applied this approximation to derive the coupled 
set of differential equations for hi (t) and Ji(t), and stud- 
ied the monopole and quadrupole excitations of bose- 
fermi mixtures after linearizing these equations around 
the equilibrium points [21]. Here we only present the 
brief derivation to make the paper self-contained. 

As specified above, we substitute the scaling ansatz Eq. 
(3) into Stringari's hydrodynamic equations. By setting 
R% = Ti/bi(t), one finds, 
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The coupled differential equations for the scaling param- 
eters bi(t) can be obtained by multiplying Eq. (9) by 
i?in° (Ri) on both sides and integrating over the spatial 
coordinates. Making use of the equilibrium properties of 
the density distribution in the ground state, 
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after some straightforward algebra one finds, 
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where (R 2 ) b = (l/N b ) J d 3 Kn b (R)R 2 is the average size 
of bosons along the z-axis. The last two terms in Eq. (11), 
linear in g b f, account for the effects of boson- fermion in- 
teraction. 

Analogous procedure can also be applied to the 
fcrmionic part. We substitute the scaling ansatz Eq. (4) 
into Boltzmann-Vlasov kinetic equation to get 
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where Ri = Ti/bi(t). Using the equilibrium properties of 
the distribution function, 

aRi 1 dVi m b dRi oVi 

to replace the first term of Eq. (12) by a linear superpo- 
sition of Ri{df /dVi) and (dnl(Ri)/dRi){df Q /dV^, and 
taking the moment of RiVi (namely, 
{1/N f ) J RiVilEq.iU^Kd^V), we finally obtain: 
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where (Rj) f = ± J ^Rn^R)^ 2 . 

The coupled set of differential equations (11) and (14) 
is a generalization of Eqs. (5) and (6) in the presence 
of the boson-fermion coupling. It determines the cou- 
pled dynamics of the condensate and degenerate Fermi 
gas in the collisionless regime as far as the assumption 
of the simple scaling solution is valid. We shall only 
interested in the evolution of bi{t) and ji(t) with initial 
condition of 6,(0) = 1, 7i(0) = 1, k(0) = 0, and 7l (0) = 
(i =± or z), in case of switching off the trap suddenly at 
t = 0, i.e., u>i b j(t > 0) = 0. Our final aim is to calcu- 
late the aspect ratio of the condensate and the Fermi gas 
defined as Xb±(t)/b z (t) and X"f±(t)/j z (t), which are ac- 
tually measured in experiment. According to Eqs. (11) 



and (14), the whole process of our numerical calculations 
in this paper consists of three stages. First, one has to 
find the equilibrium ground-state densities: n b (p, z) and 
n° (p, z) at very low temperature, which approximately 
satisfy the following coupled equations in the Thomas- 
Fermi approximation [14], 

Vho (P> z ) + 9bbn° b (p, z) + g bf n° f (p, z) = p b , 
.2/3 



2m/ 



{Q7T 2 n° f (p,z)) +Vi(p 7 z) + g bf n a b (p,z) = M/ , (15) 



where p b j is the chemical potential. It is convenient to 
obtain the solutions of Eq. (15) by iterative insertion of 
one density distribution in the other equation and numer- 
ically searching for p b and p / yielding the desired number 
of particles. Then one traces the evolution of bi{t) and 
7i(t) from t to t + At by evaluating the integrations in 
Eqs. (11) and (14), which turns out to be the most time- 
consuming step in the calculations. At the final stage, 
one computes the radial-to-axial aspect ratio. 



III. RESULT AND DISCUSSION 

In this work, we have performed numerical simulations 
with the parameters chosen to reproduce the experimen- 
tal conditions [3] . The number of bosons and fermions in 
experiment (N b = 2x 10 4 and Nf = 10 4 ) are large enough 
to ensure the validity of the Thomas-Fermi approxima- 
tion, i.e., N b a bb /a b ho JL > 1 and N f > 1 [14,22,23]. Wc 

take the harmonic oscillator length a b ho ± = \/h/(m b uj± b ) 
and hcu± b as units, and introduce the quantities a = 
rrif/m b , [3 = uj±f/cj± b , and A = uj zb /uj^ b = w z f/uj±f to 
parameterize the different mass of the two components 
and anisotropy of traps. The constraint af3 2 — 1 is al- 
ways satisfied since both bosons and fermions experience 
the same trapping potential. As in experiment, we have 
a = 0.463, = 1.47, and A = 0.0757. We have also 
taken a bb = +110ao and a b f = — 330ao, where ao = 0.529 
A is the Bohr radius [3]. The most recent measurement 
suggested a new s-wave scattering length a b f = — 410ao 
for the mixture of 40 K and 87 Rb [4,24]. However, the 
effective magnitude of a b f could be indeed lower, consid- 
ering the possible effect of the gravitational sag [3] and 
exchange correlations beyond the mean-field approxima- 
tion [25] . We will return to this point at the end of the 
paper. 

Before presenting the numerical result, it is instruc- 
tive to briefly analyze the influence of the boson-fermion 
interaction on the expansion for both species. In the 
collisionless regime, the effect of the dominated mean- 
field interaction is two-fold: (i) First of all, it influences 
the profile of the density distribution in the equilibrium 
ground state. In case of attractive interaction, as shown 
in Fig. 1, both the density of the condensate and of the 
degenerate Fermi gas are remarkably enhanced within the 
central overlap region. The condensate profile narrows 
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and the central density increases moderately. The effect 
on the Fermi gas is more pronounced: within the over- 
lap region the fermionic density exhibits a high-density 
bump on top of the low-density background and the cen- 
tral density is increased by a factor of larger than two 
[10]. Without considering the mutual attraction during 
the expansion this enhancement of the density profile, 
which corresponds to a tighter confinement, will defini- 
tively lead to a faster expansion for both species if one 
turns off the trap potential suddenly. In this paper, we 
shall consider such kind of fasten-mechanics as a static 
effect of the boson-fermion interaction on the expansion. 
(ii) On the other hand, the condensate and degenerate 
Fermi gas also interact with each other during the ex- 
pansion (especially during the early stages). As a re- 
sult, with the attractive boson-fermion interaction both 
species will feel a running confinement potential gener- 
ated by the other species and therefore reduce their ex- 
pansion rate. This slow-down mechanics is determined by 
the relative dynamics between the condensate and the de- 
generate Fermi gas. The smaller the relative expansion 
velocity is, the more attraction the two species experi- 
ence. With this respect, it is referred to as a dynamical 
effect of the boson-fermion interaction on the expansion. 

Fig. 2 shows the radial-to-axial aspect ratio of the 
condensate and of the degenerate Fermi g func- 
tion of the dimensionless expansion time variable t. Let 
us firstly concentrate on the condensate. The evolution 
of the aspect ratio with a b f = — 330ao (the dashed line in 
Fig. 2a) agrees qualitatively with the experimental result 
(the solid circles). This agreement is remarkable, consid- 
ering there is no adjustable parameter in the numerical 
calculations and the simplicity of our model. 

In order to better understand the role of static and dy- 
namical effect of the mean-field boson-fermion interaction 
on the expansion of the condensate in some details, we 
expand the density distribution of the Fermi gas around 
the center by using the fact that in experiment the Fermi 
gas distributes more widely than the condensate due to 
the Fermi statistics, and rewrite Eq. (11) into a more 
physically transparent form, 



kit) 



where 



+ 



A4 



M*)IL&i(*) 7?(*)IL7i(*) Ht)XlMt) 



o, 



9b f 



m b N b (Rf) b u>. 



- I d 3 n 



<9n°(R) 



(16) 
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Mathematically, the dynamical and static effects are rep- 
resented by the last two terms in Eq. (16), respec- 
tively. Moreover, recalling the free expansion of a pure 
condensate as shown in Eqs. (5) and (7), the quan- 
tity u>ib(l + A) 1 / 2 can in fact be interpreted as an effec- 
tive trapping frequency experienced by the condensate 
in ground state. The value of A can further be roughly 



extracted from the change of the bosonic density distribu- 
tion due to the boson-fermion attraction. The estimated 
value of A w 40% and the corresponding increase of the 
trapping frequency Alo/uj « 20% [26] is two times larger 
than that used in the fitting in Ref. [3] as we mentioned 
earlier in the beginning of introduction. This discrep- 
ancy can be easily resolved if we consider the dynam- 
ical effect of the attraction on the expansion. Indeed, 
by inserting the scaling solution for elongated traps at 
first order in A in Eq. (16), b ± (r) = yjl + (1 + S) 2 t 2 
(S represents the net increase of the trapping frequency), 
7_l(t) w y/l + /3 2 r 2 , &z(t) w 1, and j z (t) « 1, one ob- 
tains 



5 = (3 



(/3 4 + 4^4(1 + A)) 1/2 - 0- 
2A 



1/2 



- 1. 



(18) 



By substituting (3 = 1.47 and A w 40% into the above 
equation, one finds 5 = 12%, in good agreement with the 
value of 10% used in Ref. [3]. In other words, on the 
expansion of the 87 Rb condensate, the static effect of the 
mean-field attraction dominates over the dynamical one. 
If measured in units of increase of the trapping frequency, 
they contribute around+20% and —8%, respectively. 

As explicitly shown in the third term in the left side 
of the Eq. (16), the dynamical effect of the mean-field 
attraction on the expansion of the condensate is closely 
related to the relative dynamics between the condensate 
and the degenerate Fermi gas, and specifically, related 
to the value of (3 = uif/u b . In Fig. 3 the evolution of 
the aspect ratio of the condensate and of the degenerate 
Fermi gas is plotted against the dimensionless expansion 
time. For both species, the aspect ratio decreases as one 
lowers the value of /?, suggesting that the dynamical ef- 
fect of the mean-field attraction becomes more and more 
important with decreasing (3. Precisely at j3 = 1, where 
b±(r) = 7_l(t), the static and dynamical effects for the 
condensate are compensated with each other. As a result, 
the aspect ratio is almost the same as that of a pure con- 
densate. Note that although the above result for (3 ^ 1 
is not realistic for the mixture of 40 K and 87 Rb in the 
LENS experiment, it might be relevant to the system of 
6 Li- 7 Li and 40 K- 41 K, in which /3«1. 

The evolution of aspect ratio of the degenerate Fermi 
gas, on the other hand, is also determined by the com- 
petition between the static and dynamical effect of the 
mean-field attraction. As shown in Figs. 2b and 3b, 
contrary to the condensate, the dynamical effect for the 
Fermi gas always dominates and the aspect ratio is less 
than that of a pure Fermi gas. 

We now turn to consider the dependence of the as- 
pect ratio on the strength of boson-fermion interaction. 
In the first experiments at LENS [3], only the absolute 
magnitude of a b f, and not its sign, was directly deter- 
mined. It is therefore interesting to study the expansion 
also in case of repulsive interaction. In Fig. 2 the as- 
pect ratio for a b f — +330a (dotted line) is compared 
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with the result for a&/ = — 330ao, showing that the be- 
havior of the latter is closer to the experiment points. 
The fact that boson-fermion interaction is indeed neg- 
ative has been subsequently confirmed by the observa- 
tion of the mixture collapse [4]. In order to understand 
which is the general behavior in passing from negative 
to positive values of the interspecies scattering length, in 
Fig. 4 we show the aspect ratio of the condensate and 
of the degenerate Fermi gas as a function of gbf / gbb at 
a fixed expansion time. For the condensate, the aspect 
ratio has a parabolic shape with the minimum located at 
gbf /gbb = 0.5. The possible reason is that in the region 
of gbf /gbb < and gbf /gbb > 1, the condensate is always 
tightly confined by cither the attractive or strong repul- 
sive boson-fermion interaction and then the aspect ratio 
is increased by the dominant static effect of mean-field 
interaction. In contrast, the aspect ratio of the Fermi 
gas increases monotonically with increasing the value of 
gbf /gbb up to gbf /gbb ~ 2. It is interesting to note that 
in principle it is possible to tune gbf by using Feshbach 
resonances in experiment. 

Up to now, we have restricted the discussion to the col- 
lisionless regime. In the opposite hydrodynamic regime, 
the collision terms dominate over the mean-field terms 
[16]. For the condensate, to ease the analysis we shall 
still neglect the term T co u in the right side of Eq. (1), 
by assuming that the dynamics of the condensate is less 
affected by the collisions with the Fermi gas since the con- 
densate always keeps itself in the hydrodynamic regime. 
For the degenerate Fermi gas, we resort to the Euler 
equation of motion [17,18], which can be deduced from 
Boltzmann-Vlasov kinetic equation under the assump- 
tion of local equilibrium. The detailed expression of the 
coupled set of differential equations for hi (t) and 7j (t) in 
this regime can be found in Ref. [21]. The predictions 
of these equations are reported in Fig. 5. The aspect 
ratio of the condensate is affected by the collisions only 
in a minor way, which might be understood from the 
assumption of the vanishing value of T co u. For the de- 
generate Fermi gas, in contrast, the aspect ratio changes 
remarkably from the collisionlcss regime to the hydro- 
dynamic one. This is consistent with the result in Ref. 
[13], where the same Euler equation of motion has been 
used to describe the expansion of a superfluid Fermi gas 
[27]. Note that the expansion of a strongly interacting 
degenerate Fermi gas has been recently demonstrated in 
Ref. [28], and has been explained qualitatively by the 
theory of Menotti et al. [13]. However, a stringent test of 
distinguishing the system from the normal hydrodynamic 
phase to the superfluid phase is lacking in the experiment 
[29]. 

IV. SUMMARY 

In conclusion, we have studied the effect of the boson- 
fermion interaction on the expansion of a boson-fermion 



mixture, within a simple scaling ansatz. We have consid- 
ered in some details the interplay of the so-called static 
and dynamic effect of the mean-field interaction on the 
expansion. The former is caused by the modified density 
profiles in the equilibrium ground state, while the lat- 
ter refers to the interaction between the condensate and 
Fermi gas during the first stage of the expansion. These 
two effects are compensated with each other. Which one 
plays the most important role depends on the detailed 
parameters of the system. 

For the mixture of 40 K and 87 Rb, the static effect is 
dominant for the expansion of the condensate and its as- 
pect ratio is inverted more rapidly than that of a pure 
condensate. This feature has been observed in the exper- 
iment, and also well reproduced by our numerical simu- 
lations. For the degenerate Fermi gas, on the other hand, 
the dynamical effect becomes important. As a result, its 
aspect ratio is less than that of a pure Fermi gas. This 
prediction needs further experimental investigation. 

At the end of this paper, several remarks are in or- 
der: (i) The above results are based on the assumption 
of the simple scaling solution. The justification of using 
such scaling ansatz on the problem of collective excita- 
tions has been discussed in Ref [21] by the authors. We 
have shown that this approximation is equivalent to the 
sum-rule approach. However, on the problem of the ex- 
pansion, its validity deserves a further study (especially 
for the degenerate Fermi gas), though our numerical re- 
sults for the condensate are in reasonable agreement with 
the experiment, (ii) In a recent preprint, the importance 
of the corrections due to exchange-correlation energy on 
the ground state of mixture 40 K and 87 Rb has been dis- 
cussed [25]. Such corrections on the expansion can also 
be investigated by adding the exchange correlations in 
the local density approximation. (Hi) For the mixture of 
40 K and 87 Rb, the recent measurements give a more ac- 
curate value of the s-wave interspecies scattering length: 
abf = — 410a [4,24]. In Fig. 6, we show the aspect ra- 
tio for such large negative scattering length without con- 
sidering the effect of different gravitational sag for 40 K 
and 87 Rb. The result with inclusion of the exchange- 
correlation term is also plotted by the dotted line. For 
the condensate the calculated aspect ratio is larger than 
the experimental result. This discrepancy is partly due 
to the possible effect of the gravitational sag [27]. (iv) 
Finally, a careful consideration of the role of collisions is 
needed. 
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Figures Captions 

FIG. 1. The equilibrium density distribution of the con- 
densate (a) and of the degenerate Fermi gas (b) along the 
radial direction for two different values of atf. abf = (full 
line) and a&/ = — 330ao (dashed line), where ao = 0.529 A is 
the Bohr radius. We have taken abb = HOoo = 5.82 nm. The 
coordinate p and the density nbj are measured in units of the 
harmonic oscillator length a b ho± and (a b ho± ) , respectively. 
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FIG. 2. The radial-to-axial aspect ratio of the condensate 
(a) and of the degenerate Fermi gas (b) as a function of the 
dimensionless expansion time variable r with a boson-fermion 
s-wave scattering length ay = (solid line), — 330ao (dashed 
line), and +330ao (dotted line). For comparison, the exper- 
imental data (solid circles) are also plotted. One unit of r 
corresponds to 0.738 ms in the experiment [3]. Note that 
the aspect ratio for the condensate with abf — — 330ao agrees 
qualitatively with the experiment, without any adjustable pa- 
rameters. Note also that in the experiment, the typical ab- 
sorption image is taken at r = 21 (or t = 15.5 ms) for the 
condensate and at r = 6.1 (or t = 4.5 ms) for the degenerate 
Fermi gas. 

FIG. 3. Aspect ratio of the condensate (a) and of the de- 
generate Fermi gas (b) at abf = — 220ao for various values of 
j3 — u)±f/u)±b- The aspect ratio for j3 — 1.0 in figure (a) is 
almost the same as that of a pure BEC. Note that the value 
of a b j- = — 220a is not realistic in experiment. However, if we 
take abf = — 330ao, we cannot find out the density distribu- 
tions of the equilibrium ground state for f3 = 1.0 or /3 = 0.85 
within Thomas-Fermi approximation. The possible reason is 
that the mixture composed of Nb = 2 x 10 4 and Nf = 10 4 
will collapse when /3 < 1. 

FIG. 4. (a) The aspect ratio of the condensate at r = 21 
as a function of gbf/gbb- (b) The aspect ratio of the Fermi 
gas at t — 6.1 as a function of gbf/gbb- In principle, it is 
possible to tune <?(,/ by using Feshbach resonances. Note that 
gbf/gbb = 1.5875a;,/ /abb for the mixture of 40 K and 87 Rb. 
In the experiment, a b f = — 330ao corresponds to gbf/gbb = 
-4.7625. 

FIG. 5. Comparison of the aspect ratio of the condensate 
(a) and of the degenerate Fermi gas (b) at abf = — 330ao in 
the collisionless and hydrodynamic regime. For comparison, 
the result for the decoupled boson-fermion mixture in the col- 
lisionless regime are also plotted by the solid lines. 

FIG. 6. Aspect ratio of the condensate (a) and of the de- 
generate Fermi gas (b) at abf = — 410ao with and without the 
inclusion of the exchange correlation term. 
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